Magnetic Domain Wall Pumping by Spin Transfer Torque 
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We show that spin transfer torque from direct spin-polarized current apphed parallel to a mag- 
netic domain wall (DW) induces DW motion in a direction independent of the current polarity. 
This unidirectional response of the DW to spin torque enables DW pumping - long-range DW dis- 
placement driven by alternating current. Our numerical simulations reveal that DW pumping can 
be resonantly amplified through excitation of internal degrees of freedom of the DW by the current. 

PACS numbers; 75.60.-d, 75.70.-i, 72.25.-b 



Dynamics of magnetic domain walls (DW) P, 0, 0] ex- 
cited by spin-polarized electrical current are a sensitive 
probe of interactions between spin currents and magne- 
tization of itinerant ferromagnets Because DW can 
be readily detected and manipulated in magnetic nanos- 
tructures [1, @, 0, S, i, El, [H , current-induced DW 
motion provides a convenient testing ground for theo- 
ries of cou pled spin -dep endent transport and magneto- 
dynamics [l3l . 14 . 15 . [l6| . Understanding current- 
induced DW motion is also important for applications 
in nonvolatile magnetic memory [l7| . 

Most studies of DW dynamics induced by spin- 
polarized current have focused on the current-in-planc 
(CIP) geometry, in which current flows in the plane of 
a ferromagnetic film or wire, perpendicular to the DW 
plane 0, 0, i, S M, E, 0, [ll, [li. These studies es- 



tablished that long-range translational DW motion is ex- 
perimentally observed at current densities well below the 
threshold current [l^, 14 1 for excitation of DW motion 
by Slonczewski spin torque (ST) [l^, and thus other 
types of current-induced torque, such as non-adiabatic 
ST, were invoked to explain the observations 
Recent theoretical 2^, 21 1 and experimental work [22.1 2 3l| 
explored another geometry for current-induced DW mo- 
tion, in which spin-polarizcd current is applied parallel 
to the DW plane. This geometry can be realized in a 
currcnt-perpendicular-to-plane (CPP) spin valve shown 
in Fig. 1(a); where spin-polarized current from a uni- 
formly magnetized flxcd layer induces DW motion in the 
free layer [23|. 

In this Letter we theoretically examine DW motion in 
the free layer of a nanowire-shaped CPP spin valve. We 
find that Slonczewski ST from direct current is sufficient 
to induce long-range DW displacement, and that non- 
adiabatic ST jisj and field-like ST 0] are not necessary 
for translational DW motion in the CPP geometry. Our 
analytical and numerical calculations show that ST in the 
CPP geometry induces DW motion in a direction inde- 
pendent of the current polarity. This type of DW motion 
is in sharp contrast to current-induced DW dynamics in 
the CIP geometry, in which the DW velocity changes sign 
upon reversal of the current flow direction. As a result, 
alternating current (ac) in the CPP geometry induces 
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FIG. 1: (color online) (a) Schematic of a CPP nanowire spin 
valve with a transverse DW in the free layer, (b) Magneti- 
zation profile of the transverse DW in a 50 nm wide, 5 nm 
thick Py nanowire given by micromagnetic simulations for the 
applied current density J =-10 MA/cm^. (c) Micromagnetic 
magnetization profile of the DW {my{z) and mz{z)) at the 
DW center {x — 25 nm) (squares) and Eq.(l2]) fit with A as 
the fitting parameter (lines). 



unidirectional DW displacement. This DW pumping by 
ac ST originates from a non-trivial dependence of ST 
on the angle between the current polarization and the 
DW magnetization. Calculations of the DW velocity as 
a function of the ac current frequency, /, reveal that DW 
pumping can be resonantly enhanced via excitation of 
the DW internal degrees of freedom by the current. 

Motion of magnetization of the free layer in the CPP 
spin valve shown in Fig. 1 (a) is described by the Landau- 
Lifshitz-Gilbert (LLC) equation with a Slonczewski spin 
torque term [H, Hi, III: 



dm 
~dt 



—jm X H, 



a dm 
j~dt 



+ ajT{A, rn){m x p)\ (1) 



where m is a unit vector in the direction of the free layer 
magnetization, p is a unit vector in the direction of the 
fixed layer magnetization, 7 is the gyromagnetic ratio, 
Hf. is the effective field consisting of Zeeman, anisotropy 
and exchange terms, a is the Gilbert damping constant, 
aj = 2'diM ^P^'^ torque parameter, J is the current 
density, P is the current polarization, d is the free layer 
thickness, Mg is the free layer saturation magnetization, 
e is electron charge, and t(A, m) = (^Ai+i)+(A-^-i)m-p ^e- 
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scribes the dependence of ST on the angle between m 
and p [ill. In Eq.©, we omit the field-Uke ST term [i| 
because it is small in metallic spin valves [3, ISl ■ 

For a thin free layer, the DW is transverse as shown in 
Fig. 1 (b) [1, [l^l ■ The spatial profile of the magnetization 
of such a DW in the coordinate system of Fig. 1(a) is 
described by 



nix = sin(6') cos(0) = sech 
niy = sm{9) sin((/)) = sech 



A 
A 



rUz = cos{0) = tanh 



cos{<f)) (2a) 
sin(0) (2b) 



A 



(2c) 



where zq is the DW center position, (j) is the out-of- 
plane angle of magnetization in the DW center and A 
is the DW width. Figures 1(b) and 1(c) show the DW 
shape during current-induced motion in a 5 nm thick, 50 
nm wide permalloy (Py) nanowire {M^ = 800 emu/cm^, 
a = 0.01 [2^ and exchange constant A = 10~^ erg/cm) 
calculated using OOMMF micromagnctic code [2§|. The 
DW shape obtained from the micromagnetic simulations 
is in good agreement with Eq.Q. The best fit of Eq.® 
to the micromagnetic data shown in Fig. 1(c) yields A = 



22 nm [27| 



Substituting Eq.([21) into Eq.(IT]), making the rigid do- 
main wall approximation (A = const = 22 nm) and in- 
tegrating over z, we obtain equations of motion for two 
DW collective coordinates: the DW center position zq 
and the DW angle Q, [3, [2^: 



Xdt ""dt ^ " Ms ^ " 
a dzQ 



Kx) sin(20) 



dt 







(3a) 
(3b) 



where T and F will be called the torque and force drive 
terms and Kx and Ky are the nanowire shape anisotropy 
energy constants {Kx = Q.12x2ttM'^, Ky = 0.88x2t:M^) 
[30| . For current polarized along the z-axis, the expres- 
sions for T and F: 



T 



pjaj 



F 



7rcos((/)) 

can be integrated analytically to give 



t(A, iinz){l - ml)dz (4a) 
T{A^mz)(mz'my)dz (4b) 



T = pja. 



2A2 In [A] 
' A2-1 



^PiajCi{A) (5a) 



F = -TQ./ "^^^/^ tan(0) = -70,7(2 (A) tan(0) (5b) 

where p = p ■ z ~ ±1 is the fixed layer magnetization 
direction. 




FIG. 2: (color online) Dependences of the torque ^1 T (solid 
line) and force (^2 ^ F (dashed line) terms on ST angular 
asymmetry parameter A. 



For angular-symmetric ST (A = 1), F=0 and Eq. © 
reduce to the DW equations of motion derived in Ref. 
[20t . In this case, ST from direct current gives rise to 
a small reversible DW displacement and does not induce 
long-range translational DW motion . Angular asym- 
metry of spin torque (A > 1) leads to a non-zero force on 
the DW, which gives rise to a non-zero terminal DW ve- 
locity in response to direct current. The specific form of 
the ST angular dependence may be different for different 
spin valve systems 3l[ but the force term, F, is generally 
non-zero for an angular-asymmetric ST. Figure 2 shows 
that both the torque Ci(A) ~ T and the force C2(A) ^ F 
terms increase with increasing ST angular asymmetry A. 

An analytical expression for the terminal DW velocity 
under the action of direct current can be derived from Eq. 
([3]) in the approximation of small DW angle {(p << 1): 



dzo , 



7M,Ci(A)C2(A) 
2aiKy-Kx) ' 



(6) 



The quadratic dependence of the DW terminal velocity 
on current in Eq. ^ directly follows from Eq. (|5b[) stating 
that the force term F is proportional to the product of a j 
and tan((/)). Because (j) ^ aj changes sign upon reversal 
of the current polarity, zq ^ F follows. 

To understand the DW dynamics beyond the small- 
angle approximation, we numerically solve Eq.([T|) using 
OOMMF [2^. We make micromagnetic simulations of 
the DW dynamics in a 50 nm wide, 10 fiui long, 5 nm 
thick Py nanowire for the current polarization P =0.30. 
In our simulations, a current step is applied to the DW 
in equilibrium, and the location of the DW center, zo, is 
recorded as a function of time. Figure 3(a) shows zo(i) 
given by micormagnetic simulations and by numerically 
solving Eq.Q for the current densities J = ±15 MA/cm-^ 
and A = 22 nm. Eq.Q and the micromagnetic solution 
predict similar DW dynamics. We note that the sign of 
the DW terminal velocity is independent of the current 
polarity in agreement with Eq.([6]), and that for J > 
the DW velocity changes sign approximately 1.5 ns after 
application of the current step. 
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FIG. 3: (color online) (a) DW position versus time calculated 
micromagnetically (dotted and dashed lines) and in the rigid 
DW approximation by solving Eq.([3l) (solid lines) for direct 
current density ±15 MA/cm^. (b) DW terminal velocity as a 
function of direct current for two values of the ST asymme- 
try parameter A. Triangles represent micromagnetic results 
whereas lines are given by Eq.®. 



Figure 3(b) displays the DW terminal velocity as a 
function of current for two values of the ST angular asym- 
metry parameter A. Since the DW terminal velocity is 
approximately quadratic in the current density, the DW 
moves in the direction set by the pinned layer magne- 
tization vector, not by the current flow direction. This 
unusual unidirectional DW dynamics in the CPP geom- 
etry enables DW pumping - persistent translational DW 
motion in response to alternating current. Figure 4(a) il- 
lustrates this translational DW motion for three frequen- 
cies of the alternating drive current with an amplitude of 
20 MA/cm2. 
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FIG. 4: (color online) (a) Micromagnetic simulation results 
for position of the DW center versus time for three frequen- 
cies of applied alternating current with an amplitude of 20 
MA/cm^. Inset is a blow-up of a 1 ns time interval showing 
that the DW motion is a combination of oscillatory and trans- 
lational dynamics, (b) Terminal velocity of ac-driven trans- 
lational DW motion, Va, as a function of the drive frequency. 
The behavior is resonant with the resonance frequency of 7.7 
GHz. Inset shows the low-frequency behavior for A=1.4 in 
greater detail. 

Figure 4(b) shows the DW terminal velocity, Ua, av- 
eraged over one period of the drive current versus the 
drive frequency calculated (i) from Eq.(l3]) and (ii) from 
micromagnetic simulations. In the low frequency regime 
(/ < 0.1 GHz), both approaches give similar results ~ Va 
decreases with increasing frequency. This is expected be- 
cause at short time scales, the instantaneous DW velocity 
changes its sign upon reversal of the current polarity as 



illustrated in Fig. 3(a), and thus the net DW displace- 
ment is expected to be smaller for ac-driven motion than 
for dc-driven motion. However, at higher drive frequen- 
cies, micromagnetic simulations show a resonant increase 
in Va (with a resonance frequency of 7.7 GHz), in sharp 
contrast to the predictions of Eq.([3]). This indicates that 
the rigid DW approximation assuming the DW shape 
given by Eq.([2]) breaks down for high frequencies of the 
ac drive. 

We analyze the high-frequency resonant motion of DW 
driven by alternating current in the DW center-of-mass 
reference frame (COMF), defined as a the reference frame 
where zq = const =0. Figure 5(a) shows the out-of-plane 
component of magnetization, vriy, in the middle of the 
wire {x = 25 nm) as a function of time and distance 
along the wire in the DW COMF. The magnetization 
driven by alternating current at a frequency / =7.5 GHz 
with amplitude 20 MA/cm^ exhibits oscillatory behavior, 
and Fourier analysis shows that the resulting dynamics 
consist of motions at three frequencies. Figure 5(b) shows 
m,y{z, f) - the Fourier transform of my{z, t) in the COMF 
at three positions along the wire length (z). The exact 
center of the DW {z = 0, x = 25nm) is dominated by a 
GHz mode, the region 50 nm from the center is domi- 
nated by oscillations at 7.5 GHz (/), and the region far 
from the DW (z = 500 nm) is dominated by 15 GHz (2/) 
oscillations. 
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FIG. 5: (color online) (a) Out-of-plane component of magne- 
tization, my{z,t), at x = 25 nm plotted in the DW COMF 
reveals spatially antisymmetric oscillations near the DW cen- 
ter, and propagating spin waves far from the DW center for 
DW driven by 20 MA/cm^, 7.5 GHz alternating current, (b) 
Fourier transform of my{z, i), my{z, /), at three different loca- 
tions along the wire length {z = 0, 50 and 500 nm). The DW 
center {z — 0) is dominated by GHz translational DW mode. 
An antisymmetric localized DW mode with a frequency of 7.5 
GHz is dominant at z— 50 nm. Dynamics far away form the 
DW center are characterized by propagating spin waves with 
a frequency of 15 GHz and a wavelength of 65 nm. 

The DW velocity resonances with frequencies near 
GHz and 7.7 GHz are due to excitation of the lowest 
(symmetric) and second-lowest (antisymmetric) bound 
spin wave modes, which are localized in the potential well 
formed by the DW The GHz mode is the usual 

translational DW mode with cj) ; 
in the rigid DW approximation 



! [34[ that also exists 
The 7.7 GHz mode is 
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a spatially antisymmetric localized mode |33l . 34 1. Ex- 
citations at 15 GHz are propagating spin waves emitted 
by the oscillating DW at twice the drive frequency [s^ ]. 
Fig. 6 illustrates the spatial profiles of the symmetric 
{rfis) and antisymmetric {jfia) modes. 

For a drive frequency close to 7.7 GHz, the DW ex- 
hibits enhanced translational velocity, Va due to excita- 
tion of the antisymmetric mode by the alternating cur- 
rent. Indeed, Eq. (j4bp shows that the force term, F, that 
gives rise to long-range translational DW motion becomes 
large for the antisymmetric mode because ruzmy becomes 
approximately even in z, making the integral over z in 
Eq. (j4bp large. The resonant enhancement of Va takes 
place because the drive frequency must be close to the 
resonance frequency for the amplitude of the antisym- 
metric mode (and thus the force term F) to be large. 
The terminal DW velocity changes sign as the drive fre- 
quency is swept through the frequency of the antisym- 
metric mode (7.7 GHz) because the phase between the 
ac ST drive, aj, and the out-of-plane component of mag- 
netization, my, of the antisymmetric mode changes by 
180° when the drive frequency is swept through the res- 
onance frequency. The 180° phase shift leads to the sign 
change of the force term, F, in Eq. (j4bp and thus to the 
sign change of the DW terminal velocity, Va- 
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FIG. 6: Numerically determined spatial profiles of the sym- 
metric {rhs, GHz) and antisymmetric {rha, 7.7 GHz) local- 
ized DW modes. 

In conclusion, we have shown that spin transfer torque 
from direct current in the GPP geometry gives rise to 
persistent translational DW motion. In this geometry, 
the terminal DW velocity is independent of the current 
polarity, which gives rise to DW pumping - long-range 
DW displacement driven by alternating spin torque. Our 



numerical simulations reveal that DW pumping is reso- 
nantly amplified when frequency of the ac spin torque 
drive is close to the frequency of an antisymmetric local- 
ized DW mode. This work was supported by NSF Grants 
DMR-0748810 and EGGS-0701458 and by the Nanoelec- 
tronics Research Initiative through the Western Institute 
of Nanoelectronics. 
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